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Introduction 


The dimensions of space structures being considered for 
future applications are on the order of several hundred meters 
to several kilometers and will require a large number of actua- 
tors and sensors for attitude and shape control. A solar power 
satellite, for instance, may require hundreds of control moment 
gyros and thrusters to damp out surface vibrations caused by 
periodic disturbances such as solar and gravity gradient torques. 
The questions which naturally arise are: (a) where the actuators 

and sensors should be placed, (b) what typos should be used, and 
(c) how many should be used. 

Placement represents a substantial degree of freedom avail- 
able to the designer and is usually not a very straightforward 
question. It is even less apparent when one considers redundancy 
in the system to allow for failures; even if the "optimal" posi- 
tion of an actuator is known, it may not be so clear where a 
backup actuator should be placed. The answer will likely depend 
on, among other things, the operating strategy — such as whether 
or not it is intended to use all available actuators at all times. 

The types of control system components to be used is normally 
decided early in the design process based on their utility, cost, 
availability, reliability and other factors. This decision will 
not be discussed further here although the effectiveness cf 
different types of sensors and actuators can be evaluated using 
the observability and controllability measures which will be 
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developed. The number of components to be used must reflect 
the trade-off of cost, weight, power, etc. vs. system perform- 
ance— and the evaluation of performance should recognize the 
likelihood of some component failures during the lifetime of 
che system. 

In this work we develop a methodology for measuring the 
performance of a system which reflects the type, number and 
placement of the actuators and sensors on the structure. The 
measures also reflect the expected loss of performance due to 
component failures. These performance measures a_e intended to 
be especially useful as guides to the choice of component number 
and placement. 

Problem Definition 

It would be most helpful to the control engineer to have 
some criterion at his disposal for placing actuators and sensors. 
Unfortunately, modern control theory does not provide any such 
measure of "controllability" and "observability." Controllabil- 
ity is simply a binary concept— either a system is v'jontrollable 
or it is not. It does not say how controllable a system is. A 
vibratory mode of a beam, for example, is not controllable by a 
force actuator placed exactly at one of the nodes, but it is 
controllable by an actuator placed just off the node. One would 
suspect that an actuator slightly farther out would have even 
more control capability, but one can only verify that the sys- 
tem will be controllable. The same conditions hold with respect 
to observability for a sensor. 
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What should a more quantitative measure of controllability 
take into account? First, it is necessary to define a control 
objective. The most likely choice is to return the system to 
some specified state (usually the origin) after an initial dis- 
turbance. Secondly, the criterion should include how much 
control effort is required to accomplish this task. Finally, 
one should somehow standardize the criterion by the magnitude 
of the initial disturbiincf* . A larger disturbance returned to 
the origin with the same amount of control as a less perturbed 
system would likely have a more favorable degree of controllabil- 
ity. It will also be necessary to normalize the initial states 
so that one unit in each direction is equally "important," since 
rarely are all states expressed in the same units or of equal 
concern. 

Many ideas for observability parallel those for controlla- 
bility if the word "state" is replaced by "state estimation 
error" (the difference between the estimate of the state and 
the true state) : (1) the objective of measurement is to re- 

duce the error covariance toward zero, (2) accomplish this using 
the measurements optimally, and (3) standardize the criterion by 
th€. magnitude of tolerable errors. 

Previous Work 

Several papers have been encountered which deal with the 
subject of controllability and observability, but only two 
(Juang and Rodriguez [1] and Likins [2]) formulate measures using 
the types of standards just outlined. Horner [3] has considered 


optimum actuator placement but does it for the specific case 
of passive damping of a free-free beam. Skelton and Hughes [4] 
define measures in terms of controllability and observability 
"norms" which apply to the individual modes of a system rather 
than to the system as a whole. Their approach is also tailored 
to "linear mechanical systems" which have a special form of 
representation as a second order matrix differential equation. 
Although that form applies to space structure dynamics r we 
prefer to define measures which have a physical interpretation 
in terms of control or estimation error characteristics for 
general linear systems. 

In order to get a perspective on the measures of control- 
lability and observability in the sections which follow, it may 
be helpful to review the two papers which develop similar con- 
cepts. Juang and Rodriguez take an approach very similar to 
the linear quadratic regulator formulation. For the LTI state 
equation, 

x(t) = Ax(t) + Bu{t) 


they define the cost function 


J 



(x*^Qx + u'^Ru) 


dt 


where Q and R are weighting matrices on the state and control, 
respectively. This is the same cost function as for the LQ 
regulator problem except that the usual additive quadratic term 
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involving the final state is not defined because an infinite 
time horizon is allowed and xCt^) converges to zero. Thus 
the integral directly penalizes state excursion from the de- 
sired final state (the origin) as well as control effort. 

Performing the minimization on J and letting ® • 

one obtains the optimal cost function » 

J° - y x^’ct^) P°x(tQ) 

where P^ is the steady state solution of the matrix Riccati 
equation 

P - -PA - A*^P + PBR“^b'^P - Q. 

Since the control effectiveness matrix B is a function of the 
actuator locations is also a fxxnction of the actuator 

positions Thus, the optimal cost is a function of both 

initial state and actuator positions. 

For a fixed initial state, the optimal cost with respect 
to actuator positions is defined as: 

where are the actuator locations giving the minimum cost. 

Now since the initial state can have several directions in state 
space, the expectation with respect to x^ is invoked: 
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J°*{€.) - min E[J°(€)1 


or 


where 


J°*(€. ) - min i Tr(P°Q°) 
b ^2 


Q° « E[x(tQ)x(t^)'^] 

The optimal placement of actuators is then defined to be the 
position vector giving the absolute minimum of the expectation 
of the cost function . 

We found several objections to this method: 

(1) The weighting of control effort versus state excursion 
is rather arbitrary. 

(2) If there is a particular direction x^ in which the system 
is not very controllable, the information is largely lost 
when the cost is averaged over different initial states. 

(3) The degree of controllability is actually an inverse measure 
since a higher cost function represents a lower degree of 
controllaJaility and actually becomes infinite when the 
system is uncontrollable. 

(4) While control use is penalized, no effort is made to 
bound it . 

Likins develops a more sophisticated technique to be used 
in the case or bounded control effort. Using the variation of 


constants formula. 
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x(t) 




x(t. 


j ^ 


tg ,T)bu(T)<jT 


and choosing t^*0 and t»T, one can define the displacement 
in state space <5 time T 



^o " 


[I -^”^(T,0)]x^ 


+ J 5 (0,t)Bu(t)dt 


Choosing x^»0, 6 reduces to 
T 

s-!$ (0,t)Bu(t)dt » -Xq 
o 


where u of the original system has been normalized so that 
l^ij ^ 1 and B redefined appropriately. 

Likins then proceeds to define a "recovery region" 'H 
as the volume of initial i.tates that can be returned to the 
origin in time T under bounded control I I ^ 1', i.e. , 


= |x(0)j 3 'i(t) , 


t6 [0,T] , u, (t) < 1 for i»l 


/m 


x(T) 


The measure of controllability is chosen to be the minimum 
distance from the origin, over all directions in initial state 
space, of the outer surface of this region. 
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P - inf II x(0) II V x(0) ^ R 

The pr’obleni now reduces to finding the minimum norm of 

S (or x^) on this surface. This is a difficult problem which 
o 

requires, in effect, the definition of optimum bounded control 
trajectories which reach the origin in the specified time 
from many different initial conditions. Likins expresses this 
problem in terms of quadraturus which must, in most cases, be 
computed numerically. One can only compute a finite nu?nber 
of these and use the smallest computed S as the controllability 
measure. (A parallelogram approximation to the recovery region, 
such as is indicated in Fig. 1, is suggested by the authors.) 

If a system were actually uncontrollable there is no guarantee 
that one would compute the trajectory for which <5 is zero. 

The overriding objection to this method is the complication 
involved in the multiple control case. An important attribute 
of the measure of controllability will be its easy computation. 
Another objection is that Likins chooses to bound control 
magnitude and does not attempt to perform any sort of minimiza- 
tion with respect to quantity of control used, citing bounded 
control magnitude as the more realistic situation. It is 
usually the case, however, that quantity of control (e.g., fuel 
in thruster, stored angular momentum in CMG) is the primary con- 
sideration, not saturation of the controller. 
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DYNAMIC MEASURE OF CONTROLLABILITY 


The measure of controllability formulated here combines 
some of the characteristics of both of these methods. Like 
Juang and Rodriguez, it involves minimizing a cost function, 
and as Likins, the final degree of controllability involves 
a measurement in some "maximized" initial state space. The 
difference is that the cost involves only the control, where 
a quadratic is chosen for convenience to approximate magnitude, 
and the initial state is maximized with respect to integrated 
control utilization rather than running the control at satura- 
tion for the duration of the control period in question. 

The degree of controllability is the result of a four step 
procedure ; 

(1) Find the minimum control energy strategy for driving the 
system from a given initial state to the origin in the 
prescribed time. ["Control energy" is defined as 




Rudt, where R is a positive definite weight- 


ing matrix.] 

(2) Find the region of initial states which can be driven 
to the origin with constrained control energy and time 
using the optimal control strategy. This region is 
bounded by an ellipsoidal surface in state space. 

(3) Scale the axes so that a unit displacement in every 
direction is equally important to control. 

(4) The degree of controllability is a linear measure of the 
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weighted "volume" of the ellipeoid in thie equicontrol 
space. 

Step 1 can be stated mathematically as follows: 

T 

1 ^ T 

min E « y j u Rudt 
o 

subject to X ■ Ax + Bu (1) 

- x(0) - Xq 

X(T) - 0 

The Hamiltonian for this problem is; 

H - ^ u"^Ru p'^(Ax + Bu) 

so that 

P » -A*^P P(o), P(T) free 

u*(t) » 'R'^B'^P(t) 


( 2 ) 

(3) 


where u*(t) is the op:imal control. 

To find P(t), combine the differential equations (1) and 
(2) into matrix form using thr» optimal control (3) ; 
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Then donotlng the state transition matrix for the augmented 
T *p T je 

State vector [x ?“] as JXt) , and making use of the 
identities ^(0)«I and where A is the new state matrix 

in (4), the costate variable is found to bet 

where ^*espective partitions of the 
state transition matrix $t) . 

Step 2: In order to carry out step 2 rf the procedure, we 

will require an expression for the optimum cost, 

1 T 

£♦ “ j I Ru*dt, as a function of the initial state. 

To this end, we seek a relation of the form 

X » VP (6) 

since P is a function of the initial state. Differentiating 
(6) , substituting (1) , and noting that the resulting equation 
set equal to zero must hold for arbitrary P, we find that 

V » AV + VA*^ - (7) 

with the boundary condition 

V(T) » 0 (8) 

to satisfy the requirement that x(T)«0 since in general P(T) 
is not zero. We choose this boundary condition for V as a 
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matter of convenience; any other terminal value which satisfies 
the r-iquirement V(T) P(T) » 0 would produce the seune result for 
the control energy. The reason for not using the usual rela- 
tion P*Wx is that in order for P(T) not to be zero, W(t) would 
have to be poorly defined at t=T. 

Corresponding to the usual cost expression 

J = I x(0)’’w(0)x(0) 


we expect the energy cost to have the inverse form 


E = I x(Q)^V(0)"^x(0) 


(«?) 


rhe validity of this expression can be verified as follows: 


Generalize the initial time to t^. 


Then 


E - 1 
E ~ j 


I 


Rudt 


( 10 ) 


and we would like to show 


E = i x(t^)^ x(t^) 


( 11 ) 


Differentiating (10) with respect to the initial time and 
substituting (3) gives 


dE 


- j P(tQ)’^BR‘'Vp(t^) 


( 12 ) 
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Substltuting (6) into expression (11) (which is to be verified) 
we have 


E - j P(tQ)'^ V(t^) P(t^) (13) 

Differentiation of this and substitution of (2) yields the same 
result as equation (12) so that the derivative of the quadratic 
expression for E in (9) is correct. 

Also, the boundary condition matches as we can see by 
letting t^— > T. Since the optimal trajectory tends toward 
the constraint x(T)*0, the control energy E(t^) tends to 0 
as t^— ♦ T and x(t^) —>>0. The property E(t^)— ► 0 as t^— ^ T 
is assured by the form of E given in (13) and the boundary 
condition on V 

lim V(t^) - V(T) - 0 (14) 

to-T 

Equation (9) defines an n-dimensional ellipsoidal surface 
in initial state space. Any point within the ellipsoid can be 
returned to the origin in time T with energy E using the optimal 
control in eq. (3) . Though the energy expression (9) is simpler 
than that appearing in (1 ) , the differential equation for V in 
(7) remains to be solved. The solution to (7) for the case of 
rigid body and vibratory modes of a spacecraft ia presented in 
the section on Applications. 
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step 3 is to scale the axes so that a unit displacement in 
every direction is equally important. But what is meant by 
" important"? 

It may first occur to the reader to scale each state by 

the magnitude of its maximum tolerable displacement, x. , 

' ^maxl 



z 


o 





X 


o 


sc that a unit displacement in every direction is equally in- 
tolerable. But this scaling is highly inappropriate for the 
following r- son. For a fixed amount of control energy and 
time, tl larger the volume of initial states encompassed by 
the quadratic surface in eq. (9) is, the better the system 
can be controlled? larger initial states can be returned to 
the origin with the same control effort and time. Increasing 
the x^ dimension of the ellipsoid, for instance, indicates a 
favorable control capability. But if x^^ is scaled by dividing 

its maximum tolerable value, xi , we observe the following 

■^max 

paradox: as Xi is made smaller, meaning that smaller values 

■^max 
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o£ can be tolerated (or is more Important in terms of 
system performance) then the scaled variable, becomes 
larger which signifies improved control capability. 

It is apparent that the appropriate scaling should make 

a more important variable transform to a smaller value in 

the new space so as to emphasize the need to control that 

variable. The problem is that controllability should not be 

related to the accuracy with which a variable is ultimately 

controlled (which is what the above scaling does) , but rather 

to the size of the excursion one would like to be able to 

achieve. Thus let xi be the minimum state excursions one 

min 

would like to be able to return to the origin in a given time 
using a prescribed control energy. Then define the transforma- 
tion 
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so that unit values of z in any direction represent controllable 
displacements of equal importance. If controlling a given 
state is deemed less important (which is useful to recognize 
since it requires less control capability) , the corresponding 
state in z-space is made larger. 
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step 4 is to measure the controllability represented by this 
ellipsoid in equicontrol space (z-space) . Consider a two- 
dimensional case in which it is as important to control an 
initial displacement in the direction twice as large as one 
in the X 2 direction. In this case the ellipsoid defined by 
equation (9) is an ellipse in x-space. Let the ellipse have 
the shape illustrated in Figure 2a. This represents the ideal 
allocation of control since we are able to control a maximum 
displacement in the x^^ direction exactly twice as large as one 
in the X 2 direction. Figure 2b illustrates that the ellipse 
becomes a circle when transformed to equicontrol space via 
equation (15) . Thus any deviation from a circle in equicontrol 
space represents a less than ideal control allocation. 

After considering a number of alternatives » the de- 
gree of controllability was chosen to be the following; 

where is the n-dimensional volume of the ellipsoid in 
equicontrol space and Vg is the volume of the largest inscribed 
sphere; n is the dimension of the state space. The first term 
on the right side of (16) is the predominant term in the con- 
trollability measure; it reflects the smallest magnitude of 
initial state in equicontrol space which can be driven to the 
origin in the specified time using the specified control energy. 
If the controls were ideally allocated, the initial condition 
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surface would be a sphere and Vg would be the controllability 

measure. The second term in (16) adds a smaller eunount to DC 

to recognize the larger region of state space from which the 

system can recover if the surface is not spherical. The addi- 

Vg 

tional volume, scaled by ^ so that the most this 

term can add, as Vg— ► oo , is Vg and so that DC is zero if there 
is any direction from which the system cannot recover at all — 
this is the case of traditional uncontrollability, and Vg«0. 

The nth root of the weighted volume is taken as the controlla- 
bility measure to make it proportional to the linear dimensions 
of the region from which the system can recover. The volume 
weighting scheme for a two-dimensional case (volumes are areas) 
is depicted in Figures 3 (a-c^ 

Once one accepts (16) as a reasonable assessment of 
the controllability of the system, what remains to be shown are 
the mechanics of computing the n-dimensional volumes Vg and V^. 
Consider the quadratic form, x A x * d, where x is a vector of 
length n, A is an nxn matrix, and d is some scalar constant. For 
the two dimensional case., this quadratic surface is an ellipse 
and the enclosed area is given by TTab, where a and b are the 


intersections of the ellipse with its principal axes. The 

intersections are \/^ and where the A's are eigenvalues 

^ Al ^ A2 

of A so that the area equalsTTd ^ three dimensions, 

the surface is an ellipsoid and the enclosed volume is 




1 


-2 


■3 
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For n-dimenaions the volume is defined by n integrations over 
the n axes (bounded by the intersections of the surface with 
the axes) and is found to be K • where K is a 

constant. Since volume for n ^ 4 has little absolute signifi- 
cance the constant K is dropped and the volume is taken to be 
simply 


V - (V ^ 


-1 


(17) 


To apply this result to the case at hand, first sub- 
stitute (15) into (9) to obtain the eauation of the ellipsoidal 
surface in equicontrol space 

(18) 




V_ is then given by (17) where X. are the eigenvalues of (DV D) 

K X O 

From (7) and (15) we observe that both D and V are symmetric 


-1 


matrices so that the product DV^D is also symmetric. The eigen- 
values of the inverse of a symmetric matrix are just the recipro- 
cals of the eigenvalues of the original matrix. Therefore, if 
V ^ denote the eigenvalues of DV^D, the ellipsoidal volume is 
also given by 


V. 




(19) 


and the spherical volume is the shortest distance to the surface, 
1/ .yr" , to the nth power, or alternatively. 


\l ymin’" 


( 20 ) 
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The degree of controllability can then be computed using (16) , 

(19) , and (20) and actually becomes zero when the system is 
uncontrollable; the ellipsoid collapses to zero in the uncon- 
trollable direction so that y ■ is zero. 

To find the least controllable direction in equicontrol 
space (the point closest to the origin ) , we note that the 
principal axes of the ellipsoid are in the same directions as 
the eigenvectors of (DV^D)”^, and the eigenvectors of (DV^D) ^ 
are the same as those of DV^D. Therefore, the point of closest 
approach is in the direction u . , where 


DV Du . * . u . 

o min min min 


( 21 ) 


To recover the direction in the original state space, simply 

multiply u„. by D~^. 

min 

One further consideration is important in defining 
the Degree of Controllability of a system; that is how the 
measure varies with number of actuators. The Degree of Con- 
trollability has been defined in terms of a constraint on 
control energy with no reference to a constraint on control 
magnitude. But it seems appropriate to recognize the fact 
that a system with more actuators has greater control capability 
when there is a limit on control magnitude — as is always the 
case. The measure of controllability as defined above can be 
made to vary directly with the mamber of actuators placed at 
the S 2 une locations by scaling the elements of R inversely with 
m — ^the number of actuators in the sys+‘'“m. Usually R is taken 
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diagonal , and if the diagonal elements are first chosen to 

reflect the relative cost of using the different actuators, then 
the final elements of R are defined to be 

Rii - Rq /m (22) 

with m « total number of actuators. 

Dynamic Measure of Observability 

Any measure of the observability of a dynamic system 
should reflect as directly as possible the amount of information 
which can be derived about the system states from the sensor 
outputs in a given aunount of time. The means of obtaining this 
information is by attaching to the system an observer whose 
states, X, are "estimates" of the true states of the system. 

The more information that is obtained about the system, the 
smaller the estimation error becomes. 

A direct indicator of the amount of information one 
has about the system states is the information matrix, the in- 
verse of the error covariance matrix. In order to maximize 
the amount of information, one should minimize the estimation 
error. The linear estimator which minimizes the state estima- 
tion error vector, e = x - x, in a mean square sense, i.e., 
minimizes 

S - e*^ Me (23) 

where M is some weighting matrix, is the Kalman Filter. 

For the Kalman Filter, the error covariance equation 




is 
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P * AP + PA*^ - PC^n"^CP + Q (24) 

where P is the estimation error covariance matrix, and JI and Q 
are the measurement and driving noise intensity matrices, respec- 
tively. Since the measurement noise is a property of the set of 
sensors being evaluated, we retain its inclusion in (24) in the 
form of N but do not include the effect of state driving noise, 
because that is an external influence not related to the sensor 
set. Thus, if we set Q»0 and call the information matrix 
J(*P'"^), then (24) in terms of J becomes 

J » - JA - a'^J + (25) 

Take as the standard situation the case in which there 
is no information about the state initially and data is collected 
up to a specified time T. Then J(0) -*0 and one is interested 
in J (T) . Having the information matrix at time T, we are inte- 
rested in measuring how much information has been accumulated. 

One way of measuring the size of J(T) is by reference to the 
quadratic surface 

j“^v . 1 (26) 

As with equation (9) in the control case, equation (26) defines 
an ellipsoidal surface in v-space. If J is a diagonal matrix 
(one can always transform to principal coordinates) , one observes 
that increasing an element will expand the ellipsoid in the 
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direction Thus the larger J becomes, the larger the 

volume encompassed by the surface in (26) so that the more 
information obtained about the system, the larger the volume 
becomes . 

Typically, however, some components of x will be 
of greater concern than others — especially considering that 
different units will apply to different components. Paralleling 
the discussion of the control case, define the transformation 

w ■ Fv 


F 





1 

max 


(27) 


where e. are the maximum errors one is willing to tolerate 
^max 

in the direction x^. The more error one is willing to tolerate 
in that direction, the greater the transformed state so the 
larger the volume becomes. Thus the scaling is consistent with 
the ideas presented in the last section. Also note that v has 
units of reciprocal error, so w is dimensionless as was z in 
the control case. 

Now that the axes have been scaled so that it is 
equally important to obtain information in each direction, one 
can use the same definition for the degree of observability as 
was used for controllability when applied to equicontrol space. 



- 23 - 


Agaln, the ideal aensor distribution would produce a sphere 
in 'v-space^ so that the degree of observability involves a 
spherical volume plus a lesser weighted excess volume due to the 
nonideality of the distribution. Specifically, 


DO 



(28) 


with 



and the V ^ are the eigenvalues of FJ(T)r. 

The remaining problem is to solve the differential 
equation (25; for J so as to write out explicitly J (T) . We have 

* T T -1 

J = - JA - A^J + C^N -^C 


J(0) - 0 


This is similar to the corresponding problem in the definition 
of the degree of controllability. There we required V(0) with 

• <P -IT 

V - AV VA‘ - BR 


V(T) - 0 

Define a backward time variable, T« T 


t. 


so that 


dJ 

37 


dJ 
■ It' 


Then in terms of T» equation (25) becomes 

• T T -1 

J - JA + A^J - C N 

J(T) -0 
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(29) 


This is the same as the equation and boundary condition for V 
wi\,h the substitutions: 



R N 


So if a subroutine is prepared to produce V(0) given A, B, and 
R, that same subroutine can be used to produce J (T) by use 
of the substitutions indicated. 

It is worthy to note that the parallelism in computing 
the degrees of controllability and observability stems from the 
similarity between the quadratic forms (9) and (26), respectively. 
However, the concepts which drove us to those forms were quite 
different. Equation (9) represents an actual ellipsoid in 
state-space which bounds the initial states that can bo returned 
to the origin in time T with a prescribed energy £. For the 
observability case, the information retrieval capability is 
alre^idy maximized through the use of a Kalman Filter, and one 
is simply trying to formulate a measure of observability based 
upon the size of the final information matrix. Thus equation 
(26) serves only as an aid to the definition of the size of J, 
and the space in which it is defined serves only to measure that 
size volumetricaily . 
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APPLICATION TO ONE -DIMENSIONAL CASE 

To dORionfttrate the procedure for obtaining the degree 
cf controllability and observability, the above results were 
applied to the vibratory inodes of a free-free beam. Start with 
a series expans ‘on for the beam displacement y, 

y(€,t) - 

i 

where ) is an orthogonal set of modal shapes and are 

the modal eunplitudes , and substitute this into the governing 
differential equation for a beam 


El 



m 



f (6 .t) 


where f is the forcing term and m, E, and I are the beam mass 
(M) /length (£) , modulus, and cross-section inertia, respectively. 
Assuming the use of m point force actuators. 


m 

f ( € ,t) - 16 ( € - Cj) Uj (t) 

j-i 


with € j being the actuator positions and u^ (t) the control 
magnitudes, one obtains the relations 


j-1 


( 30 ) 
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where(*/j^ is the frequency of the ith mode. 

The modal shapes for a free-free beam are given by 

^l(x) = 1 



(31) 


^^(x) = cosh^^x + cos ^x - a^(sinh^j^x + sin /?ix) i>3 

where the are the solutions to 
1 - cosh cos = 0 

and 

sinh + sin J3^Z> 

^ coshj2 - cos P^Z, 

The first two modes of the beam are rigid body modes and thus 
have a frequency equal to zero. has the interpretation of 

the rigid body translation of the center of mass of the beam, 
and represents rotation of the beam about its center of mass. 
Next, consider casting (30) into the state space 

form, 

X = Ax + Bu 

( 32 ) 


y = Cx 
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where 
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where the number of modes has been truncated at N, and the use 
of M force actuators at positions €j and P translation rate 

at positions tti has been assumed. The replacement of 


sensors 
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a force actuator at by a torque actuator would involve 
replacing the corresponding elements of B by 
for i » 1,...,N. The use of a deflection sensor at would 
involve switching 0 and in each of the pairs 

[0 0j(CX^) ] in the ith row of C. To include natural damping 
in the model, the negative of the damping term, 2 would 

appear in each diagonal block of the system matrix of (32) multi- 
plying theXj^ term. For the present, this is considered negligible. 

Equation (V^ remains to be solved before the degrees 
of controllability and observability can be computed. The solu- 
tion of this equation is facilitated by use of the following 
real invertible transformation: 



T 




(33) 


where the v^ are the generalized eigenvectors corresponding to 
the zero eigenvalues and the a^, b^ are the real and imaginary 
parts of the eigenvector corresponding to the complex eigenvalue 
A.i - j 

If a new matrix M is defined by the relation 
V « TMT*^ (34) 

and .A is formed from the eigenvalues , 



- 29 - 



then substitution of both of these relations into ( 7 ) yields 

M = AM + MA*^ - (36) 

This equation is much simpler to solve than equation (7) for V, 
and the solution for M is presented in Appendix A. Conversion 
back to V is attained through use of (34) . 

A computer program was written to calculate the degree 
of controllability (observability) for up to four actuators 
(sensors) placed ar various positions along a free-free or 
simply supported beam (FORTRAN listing appears in Appendix B) . 

The programmer specifies the number of equally spaced positions 
along a half beam length to be tested (mode shapes are symmetric) , 
and the program computes the degree of controllability for all 
possible arrangements of actuators. The same program is used to 
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compute observability with the appropriate changes outlined 
in the last section. The present program assximes the use of 
force actuators or translation rate sensors but can be easily 
modified for torque actuators and deflection sensors. 

The program accepts as input the system matrix A, 
the number of flexible modes to be considered (maximum 5) , 
the number of actuators to be tested, the input weighting and 
control scaling matrices R and D, and the control period T. 

The mass, length, and modal frequencies of the beam were 
chosen to correspond to those of the experimental beam set up 
at NASA Langley Research Center (£i = 12 ft, m = 0.50 slugs, 

* 11.47 rad/sec, = 31.63 rad/sec.) In all trials, there 
was no relative weighting of actuators (R = I) , and the 
amplitude rates were scaled by 1/{J ^ relative to their respec- 
tive amplitudes using D (eunplitudes were considered equally 
important) . 

In Figures 4 and 5 , the degree of controllability 
(DC) is plotted for one force actuator varied along the length 
of a single mode beeun. Figure 4 shows the expected correspondence 
between the DC and the first mode shape. The maximum DC is 
at the ends where there is maximum deflection, and the DC be- 
comes zero at the nodes where the system is uncontrollable. 

The correspondence between mode shape and degree of controllabil- 
ity is again apparent in Fig. 5 when the second mode is con- 
sidered alone. 



- 31 - 


Figures 6-8 consider the first and second inodes 
simultaneously. In Fig. 6, a single actuator is tested 
along the length of the beam as in the previous two cases. 

The maximum DC is again at the ends but the system becomes 
uncontrollable at a node of either mode. The DC has an 
intermediate peak at the 7th test position which corresponds 
to an antinode of the 2nd mode. 

In Fig. 7 one actuator is fixed at the middle of 
the beeun (antinode of 1st mode) while the other is varied. 

There is an overall increase in controllability because 
of the presence of the second actuator, but the DC still goes 
to zero at the nodes of the second mode because the fixed 
actuator is at a node of the 2nd mode and thus contributes 
nothing to the controllability of that mode. The degree 
of controllability never goes to zero in Fig. 8 when the 
fixed actuator is at the end. The optimal placement of the other 
was found to be at position #7 if duplicate positioning at #1 
is not allowed. 

The degree of observability (DO) for two cases is 
illustrated in Figures 9 and 10. In Figure 9, a rate sensor 
was varied along the length of a single mode beam. The re- 
sultant DO is strikingly similar to the DC of Fig. 4. The 
first and second modes are considered in Fig. 10 where one 
sensor is fixed at the center of the beam and the other is 
varied. The DO becomes zero at three points because the second 
mode is unobservable at the location of the first sensor. 
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CONCLUSION 

While it is difficult to consider the degrees of 
controllability and observability just developed in an abso- 
lute sense, they serve well as quick relative measures of 
controllability amd observability. A more realistic measure 
of controllability, for instance, might involve the integral 
magnitude of control effort rather than the integral quadratic 
form chosen for convenience. This degree of realism has been 
sacrificed in favor of the analytic solution to the optimal 
control problem. It is also true that the "size" of the 
information matrix could have been defined in several other 
ways, e.g. , tr J, in computing the degree of obseirvability . 

The control period is also somewhat arbitrary, but if the modal 
periods are short compared to T, the measures of controllability 
and obseirvability are independent of T in a relative sense. 

The control measure does have several advantages over 
the methods in [1] and [2] : (a) it does not arbitrarily weight 

state excursions against control effort, (b) it calls attention 
to the most uncontrollable direction by primarily weighting the 
volume generated by that minimum distance — thus it is a worst 
case analysis, (c) it seeks a control law minimizing integrated 
control use, and (d) it is relatively s.v'^'rle to compute. 

For the observability case, the Kalman Filter already 
provided the minimized least square estimate error for which 
the covariance matrix is P. P determined the information matrix 
J whose size was used to compute the degree of observability. 
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The choice of measuring the size or J by the weighted volume 
within a quadratic surface made the computation of observa- 
bility analogous to controllability. 

The results of the DC and DO calculations in the 
case of the free-free beam were entirely intuitive and could 
have been anticipated from knowledge of the mode shapes. But 
that exeunple was taken in order that one could interpret the 
results easily. The purpose in defining these measures of 
controllability and observability is to assist the designer 
of a control system for a plant of realistic complexity where 
the best locations of sensors and actuators may not be so 
obvious , 

Now that these tools have been developed, they will 
be applied to the problem of choosing the number and location 
of sensors and actuators in the design of a large space 
structure considering the likelihood of random failures among 
these components. It is expected that the optimum locations 
for components with possibility of failure will differ under 
certain circumstances from those with no chance of failure. 
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A PPENDIX A 

SOLUTION OF THE MATRIX DIFFERENTIAL EQUATION (36) 

This Appendix presents the solution to the differential 
equation 

M » -A.M + - D 

where A is given by (35) and the driving matrix D is the 
last term in (36) . 

The solution matrix M(t) is symmetric and has the follow 
ing form: 


^^W = 


II 

II 

II 

... II 

II 

II 

II 

... II 

III 

IV 

IV 

... IV 


III 

IV 

... IV 



III 

... IV 




III 


The Roman Numerals indicate 2x2 block solution types. If 
the two rigid body inodes are not included in the model, the 
first and second row and column blocks are deleted from (A-2) 


(A-1) 


(A-2) 


The block solutions have the form 


If the solution 
Note that a and 


A- 2 


rruc mc4 
mb4 

is symmetric «» "'ad^ ' only m^^ is given, 

b are row indices, c and d are column indices. 


A- 3 
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Appendix B 


FItuC; OCOCON FORTRAN A 


VM/SF CQNVCRSATir^At MONITOR SYSTCM 


THIS PROGRAM COMPUTES THE DEGREE OF 
CONTROLLABILITY ANO OBSERVABILITY 
FOR FORCE ACTUATORS ANO RATE SENSORS ON A BEAM 


INPUT: N - NUMBER OF SYSTEM STATES 

NA • NUMBER OF ACTUATORS (SENSORS) 

IPHI - (1) FOR FREE-FREE BEAM 

(3) FOR SIMPLY SUPPORTED BEAM 
rSTART - FIRST ACTUATOR TEST POSITION 
<«POS - NUMBER OF POSITIONS TO BE TESTED 
(AS - (1) TO COMPUTE CONTROLLABILITY 
(2) TO COMPUTE CBSERVABILITY 
IFIX - FIXED POSITION OF SECOND ACTUATOR WHEN 

PLOTTING CONTROLLABILITY FOR 3 ACTUATORS 
8M • BEAM MASS 
BL - BEAM LENGTH 
OT - CONTROL PERIOD 

FB - FRACTION OF BEAM LENGTH FROM END OVER 
WHICH ACTUATORS PLACED 
TOL - ZERO TOLERANCE FOR REAL NUMBERS 
TAU - ACTUATOR MEAN TIME TO FAILURE 
TOP - SYSTEM OPERATING OR MISSION PERIOD 
OM - BEAM MODAL FREQUENCIES 
BETA - MODAL SHAPE PARAMETERS 

DOIAG - DIAGONAL ELEMENTS OF STATE WEIGHTING MATRIX 
R - ACTUATOR WEIGHTING MATRIX 
A - SYSTEM MATRIX 

OUTPUT LOC ' EIGHT DIGIT LOCATION CODE REPRESENTING 
POSITIONS OF 4 ACTUATORS: RIGHTMOST 

PAIR REPRESENTS LOCATION OF FIRST ACTUATOR 
ANO leftmost the FOURTH ACTUATOR (IF THE PAIR 
EQUALS IFOOOE-NOPOS^I. THE ACTUATOR HAS FAILED 
LMAX - ACTUATOR LOCATIONS FOR MAXIMUM OC NOT 
CONSIDERING FAILURES 

OCMAX - MAXIMUM DC NOT CONSIDERING FAILURES 
LMAXF - ACTUATOR LOCATIONS FOR MAXIMUM AVERAGE OC 
OCMAXF - MAXIMUM AVERAGE OC (FAILURES CONSIDERED) 
UMIN - LEAST CONTROLLABLE DIRECTION IN ORIGINAL 
STATE SPACE ASSOCIATED WITH MAXIMUM OC 
UMAX - MOST CONTROLLABLE DIRECTION IN ORIGINAL 
STATE SPACE ASSOCIATED WITH MAXIMUM OC 


DIMENSION A( 10. 10). S( 10. 4 ) . R( 4 .4 ) . I ACT(4 ) ,0M(5 ) . 

« V( 10. 10). C 1(24). WK 1(55.9). AA( 10, 10), 

i o( 10 . io).oo:ag( io).ov( to. io).ovo( io. to).cv( io). 

A WK(2OO).0VOSVM(5S).UMlN(1O).UMAX( 10). 

A RACT(4).WKAREA( 10) . RINV(4 , 4) ,BRINVB( 10, 10). 


•0EG00010 

0EQOOO3O 

0EG00030 

0EGO0O4O 

OEG00090 

OEGOOOeo 

OEG00070 

OEQOOOaO 

0EQ0(X>9O 

OEQOOtOO 

0EQOO11O 

DEQ00120 

0EQOO13O 

0EG00140 

0EQ0O1SO 

DEG001B0 

0EG00170 

DEG00180 

DEG00190 

0E000200 

0EG00210 

DEG00220 

OEG00330 

0EGOO24O 

DE0002SO 

0EC00260 

0EGOO37O 

DEG00280 

0EG0029O 

0EGO03OO 

OEG00310 

0EGOO32O 

0EGO033O 

OEQ00340 

DEG00350 

DEG003B0 

OEG00370 

OEG00380 

0EG00390 

OEQ00400 

DEQ00410 

OEG00420 

0EG00430 

OEQ00440 

DEG00450 

0EG00460 

0EGOO47O 

0EGOO48O 

DEC00490 

OEGOOSOO 

0EGOOS1O 

0EG00530 

0EGOO53O 

0EOOO54O 

OEGOOSSO 


o o u o c, o 
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PtLE: OEGCQN FORTRAN A 


VM/SP CONVERSATIONAL MONITOR SYSTEM 


RINV8T(4. tO).OC4( 12. 12. 12. 12),0C3( 12. 12. 12).DC2( 12. 
DCK 12).8ETA(5).Z( 10. 10).T(10. 10).TINV(10. 10). 
TTRANI 10, 10).TTINV( 10. 10).8RI8TT( 10. 10). 

0R( 10. 10).RW2(20).RZ2(200).AM( 10. 10),AMTT( 10. 10). 
DVOINVl 10. 10).WKAR2(5S).XARRAV(23).YARRAV(23) 
COMPLEX W2( 10).Z2( 10. 10).ZN 
EQUIVALENCE! W2( 1),RW2( t)).(Z2(1. 1).RZ3( D) 

OATA IN.IO.IOCT. tNO.NW.IJOB.EPS/5.6.0. 1.SS. 1. 1.E-1S/ 

CALL PL0TS(IDUM.IDUM.9) 

• READ AND ECHO INPUT 

REA0(IN.4) N.NA. IPHI .ISTART.NOPOS. IAS. IFIX 
REAOdN.5) 8M.BL.0T.FB.T0L.TAU.T0P 

READ! IN. 6) (OM(I).I*1.Si.(BETA(I),I-1.5).(OOIA6(I).I*i. 10) 
REAOdN.7) ((R(I.J).J-1.4).I«1.4) 

REA0dN.8.EN0-17) ( (4(1,0). J-1 ; T-1,10) 

4 F0RMAT(7I2) 

5 F0RMATOF10.4/4F10.4) 

6 F0RMAT(3(5F10.4/).SF10.4) 

7 F0RMAT(3(4F1O.4/).4f 10,4) 

8 FORMAT! 19(5F15.4/).SP1S.4) 

17 WRITE(I0,20) N.NA. IPHl. IAS, IFIX, ISTART.NOPOS. BM.8L.DT. 
a F8,TAU.T0P.(0M(I).I-1.5).(BETA(I).I-1.5). 

a ((R(I.0>.J>1.4).1«1.4).(0D1AG(I).I-1. 10) 

20 FORMAT! IX. 'N- ' . I2/*NA»* , 12/* IPHI«» . 12/' 1AS-' . 12/' IFIX-' . 12/ 
a 'ISTART-'.12/*N0P0S»'.12/ 

a '8M-',F10.4/'8L-* .FtO.4/'0T-*,F1O.4/'FB-' .F10.4/ 

a 'TAU-'.E15.4/*TOP-'.E15.4/'OM!1-S)-'.SF10.4/ 

a 'BETA! 1-9)-'. SF10.4/ 

a 'R-'/4!4F10.4/)//'00IAG(1-1O)-'.5F1O.4/12X.5F10.4////) 

IFCOOE-NOPOS+1 

• INITIALIZE VARIABLES 

00 23 1-1. S5 
DV0SYMd)-0. 

23 CONTINUE 

DO 24 1-1. 12 
0C1(I)-0. 

DO 24 0-< . 12 
OC2!I.O)-0. 

00 24 K- 1.12 

OC3d.O.K)-0. 

DO 24 L-1. 12 
OC4d.O,IC.L)-0. 

24 CONTINUE 

00 29 1-1.10 
00 29 0-1.10 

ovd.o)-o. 

Zd.0)-0. 

ovod.o)-o. 

BR1NVB!I.O)-0. 

29 CONTINUE 

00 36 1-1.10 


12), OEG00560 
0EGO057O 
0EG00980 
0EGO059O 
OEG00600 
0EGO061O 
0EGO06LO 
0EGO063O 
OEG00640 
0EG006S0 
0EGO066O 
OEG00670 
0EG0068O 
0EG00690 
DEG00700 
0EG00710 
0EGOO72O 
0EGOO730 
OE600740 
0EGOO75O 
OEG00760 
0EG00770 
0EGO078O 
OEG00790 
0EG00800 
0EG00810 
0EGOO82O 
0EGOO83O 
0EGOO84O 
0EGOO8SO 
0EGOO66O 
0EGOO87O 
0EGO088O 
0EG0O89O 
0EGOO9OO 
0E6OO91O 
0EGOO92O 
0EGOO93O 
0EG00940 
DEG009S0 
0EG00960 
OEG00970 
DEQ00980 
0EG0O99O 
OEG01000 
0EQO101O 
0EG0102O 
0EGO1O3O 
0EG01O4O 
0EG010S0 
0EG0106O 
OEG01070 
0EG01O8O 
0EGO1O9O 
0EQ01100 


o o o r> net n 
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FILE: OEGCON FORTRAN A VM/SP CONVERSATIONAL MONITOR SYSTEM 



0(I.I)-OOIAG(n 

0EGO1 1 10 


00 36 U- 1 . 10 

0EQO112O 


IF(I.NE.O) 0(I,U)-0. 

0ECK3113O 

36 

CONTINUE 

DEG01140 


NE-N*(N>1)/2 

DEG011S0 


TENO-OT 

OEG01160 


0CMAX«0. 

OEG01 170 


OCMAXF-O. 

0EQO1 180 


LMAX-O 

OEG01190 


LMAXF-0 

0EGO13OO 

0EGO121O 

FIND THE TRANSFORMATION MATRIX T USED IN COMPUTING V 

0EGO122O 



0EGO123O 


DO 671 I>1,N 

0EGO13<.v 


DO 671 J>1.N 

0EQ012S0 


AA(I.J>*A(I. J) 

0EGO1360 

671 

CONTINUE 

0EQO127O 


CALL EIGRF(AA.N. 10. I JOB. RW2. RZ2. 10.WK.IER) 

DEG01280 


DO 672 I*1.N 

0EGO129O 


DO 673 J-1.N.2 

OEG01300 


T(1.0)-REAL(Z2(I,J)) 

0EG0131O 

672 

CONTINUE 

OEG01320 


DO 673 I-1.N 

0EQ0133O 


DO 673 J>3,N.3 

0EG0134O 


T(I.J)-AIMAG(Z2(I.J-D) 

0EG013SO 

673 

CONTINUE 

0EQ0136O 


CALL LINV1F(T.N. 10. TINV. lOGT .VIK. lER ) 

DEG01370 


DO 674 I-1.N 

0EGO136O 


00 674 J-1.N 

DEG01390 


TTHAN(I.J)-T(J.I) 

DEG01400 

674 

CONTINUE 

DEG01410 


CALL LINV1F(TTRAN.N. 10. TTINV , IDGT . WK. lER ) 

DEG01420 

DEG01430 

FOURTH ORDER OO-LOOP TO PERMUTE LOCATIONS OF 4 ACTUATORS 

DEG01440 

(NO TWO locations ARE ALLOWED TO BE THE SAME) 

DEG014S0 



DEG01460 


00 46 1-1.4 

0EG01470 


lACTd )-IFCOOE 

DE001480 

46 

CONTINUE 

DEG0149O 


1ACT4-IACT(4) 

DEG01500 


IACT3-IACTO) 

0EQ01S10 


IACT2-IACT<2) 

DEG01S20 


IACT1-IACT( 1) 

0EGO153O 


IF(NA.NE.4) GO TO 49 

0EG01S40 


DO 161 :ACT4-1START.IFC0DE 

DEQ01SS0 


IACT(4)-IACT4 

D -GO 1560 


00 TO SO 

D-G01570 

49 

IF(NA.NE.3) GO TO St 

DEG01S80 

SO 

DO 171 IACT3-ISTART, IFCOOE 

DEQ01S90 


IACT(3)-IACT3 

0EGO16OO 


GO TO S3 

0EOO161O 

SI 

IF(NA.NE.2) GO TO S3 

DEG01620 

S3 

00 161 IACT3-ISTART. IFCOOE 

0EQO163O 


IACT(2)-IACT2 

0EG01640 

S3 

DO 181 lACTI-ISTaRT. IFCOOE 

0E0O16SO 


onnnno non non 


PILE: OCGCON FORTRAN A 


VM/SP CONVERSATIONAL MONITOR SYSTEM 



IACT(1)-IACT1 

DE60166O 

DE601670 

> COMPUTE CONTROL EFFECTIVENESS MATRIX B 

0E60168O 

OE601690 


DO 62 1-1.10 

0E6O17OO 


00 62 J-1.4 

0EGO171O 


B(I.J)-0. 

DE601720 

62 

CONTINUE 

DE601730 


00 63 1-2. N, 2 

0EG01740 


00 63 J-1.NA 

0E601750 


RACT ( J ) ■ ( FLOAT ( I ACT ( J ) - 1 ) / F LOA T ( NOPOS - 1 ) ) -BL • FB 

DE601760 


IF(IPHI.E0.2) 60 TO 625 

OE601770 


B(I.J)-PHI(RACT(J).BET;.(I/2).BL)/BM 

0EG01780 


60 TO 627 

DE601790 

625 

B(I.J)-PHI2(RACT(J).I/2.BM.BL)/BM 

DEG01800 

627 

IF(IAS.E0.2) B(I.J)-B(I.J)-BM 

DEG01B10 

63 

CONTINUE 

0EG01820 

DE601830 

• ZERO-OUT COLUMNS OF B ASSOCIATED WITH INOPERATIVE ACTUATORS 

OEG01S40 

DE601850 


IF(IACT4.NE.IFC00E) 60 TO 633 

DE601860 


00 632 1-2. N. 2 

0E6O1870 


8{I.4)-0. 

0E6O188O 

632 

CONTINUE 

0 EGO 1890 

633 

IF(IACT3.NE.IFC00E) 60 TO 635 

0EGO19OO 


00 634 1-2. N. 2 

0EG01910 


B(I.3)-0. 

0EGO192O 

634 

CONTINUE 

0EGO193O 

635 

IF(IACT2.NE.IFC0DE) 60 TO 637 

DEG01940 


DO 636 1-2. N. 2 

DEG01950 


B(I.2)-0. 

0EGO196O 

636 

CONTINUE 

OE601970 

637 

IF(IACTI.NE.IFCQOE) 60 TO 65 

DEG01980 


00 638 1-2. N. 2 

OEG01990 


8(1. 1)-0. 

DEG02000 

638 

CONTINUE 

OEG02010 

65 

N8-0 

DE602020 


00 66 1-1, N 

0EGO2O3O 


00 66 J-1.NA 

0EGO2O4O 


IF(ABS(B(I,J)).LT.TOL) NB-NB-M 

0EGO2O5O 

66 

CONTINUE 

DE6O206O 


IF(NB.EO.N-NA) 60 TO 151 

0t602070 

OEG02080 

' IF ALL 

ACTUATORS INOPERATIVE. 60 TO NEXT TEST LOCATION 

0EGO2O9O 

OE602100 


I TOTF - I FCOOE • 10* -B* I FCOOE • 10* • 4* 1 FCOOE • 100+1 FCOOE 

OEG021 10 


LOC-I ACT4*10*-6+l ACTS* 10**4+I ACT2* 100+1 ACT 1 

0EGO212O 


IF(LOC.EO.ITOTF) 60 TO 203 

0EGO2 1 30 
0EGO214O 

ADJUST 

INITIAL R TO ACCOUNT FOR ACTUATOR SATURATION 

0E6O215O 

0EGO216O 


NOA-O 

0EGO217O 


DO 661 1-1.4 

0EGO218O 


IF(IACT(1).NE.1FC00E) NOA-NOA+1 

DEG02190 

661 

CONTINUE 

0E6O22OO 
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FILE: OEGCON 

FORTRAN A VM/SP CONVERSATIONAL MONITOR 

SYSTEM 



00 663 1-1.4 

0EGO221O 



RINVd. I)-FLOAT(NOA)/R(I.I) 

0EGO222O 



00 663 J-1.4 

DEG02230 



IF(I.NE.J) RINVd. J)-0. 

0EOO224O 

663 


CONTINUE 

0EGO32SO 

C 



0EG02260 

C« COMPUTE 

ORIVINQ MATRIX IN O.E. FOR M 

0EQ0227O 

C 



0EGO228O 



CALL VMULFP(RINV.B.NA.NA.N 4. tO. RINVB . 4 , I ER ) 

0EQO229O 



CALL VMULFF(B. RINVB. N.NA.N. 10.4.BRINVB. 10. lER) 

0EGO23OO 



CALI VMULFF(BRINVB.TTINV,N.N.N. 10. 10.BRIBTT. 10.IER) 

0EGO231O 



CALL VMULFF(TINV.BRI8TT,N.N.N. 10. 10, OR. 10.IER) 

0EGO232O 

c 



DEQ02330 

C* COMPUTE 

DIAGONAL BLOCKS OF M (TYPE III) 

0EGO234O 

C 



0EGO235O 



"0 675 1-1, N. 2 

0EQ0236O 



SIQ1-REAL(W2(I)) 

0EQ02370 



0M1-AIMAQ(W2d)) 

OEG03380 



CALL 0IAG(0T,SIG1.0M1.0R(I. I ).0R(I.l4'1).0RdT 1.14^1), 

0E00239O 

6 


AM ( I . I ) , AM d . 1 4 1 ) . AM ( I <•> 1 . 1 4- 1 ) ) 

0E0024OO 



AM(l4l.I)-AM(I.l4t) 

0EQ02410 

675 


CONTINUE 

0EQ0242O 

C 



0EG0243O 

C* COMPUTE 

OFF-DIAGONAL BLOCKS OF M (TYPE IV) 

0^002440 

C 


, 

0EGO24SO 



IF(N.LT,4) GO TO 70 

0EG0246O 



NM3-N-3 

0EG0247O 



NMt-N-1 

0EG0248O 



00 676 I-1.NM3.2 

0EQ02490 



IP2-I43 

OEQ02500 



00 676 J-IP2,NM1,2 

0EG0251O 



SIG1-REAL(W2(D) 

0EGO252O 



0M1-AIMAQ(W2(I)) 

0EQ0253O 



SIG2-REAL(W2(U)) 

0EQ02540 



0M2-AIMAG(W2(U)) 

0EGO2S5O 



CALL 0FDIAG(0T,SIGt.SIQ2.0M1.0M2.0R(I.J). 

0E002S6O 

a 


0R(I.j4l).0R(l41.J).0Rd41.14l). 

OEG02570 

& 


AM(I.J).AM(I. J4-1).AM(I4;.J). 

0EQ02580 

& 


AM(I41,J41)) 

0EQ0259O 



AM(0,I)-AM(I.J) 

DEQ02600 



AM(U41.I)-AM(I,J41) 

OE002610 



AM(d.I41)-AM(l4l.O) 

0EG0262O 



AM(U4-1 ,I4>1)-AM(I41,J41) 

DEQ0263O 

676 


CONTINUE 

0EQO264O 

C 



0E00365O 

C* TRANSFORM 

FROM M TO V 

0E002660 

C 



0EQ0267O 

70 


CALL VMULFF(AM.TTRAN.N.N.N. 10. 10.AMTT. 10.IER) 

0EQ0268O 



CALL VMULFF(T.AMTT.N.N.N. 10. 10. V. 10.1ER) 

DEG02690 

C 



DE002700 

C* TRANSFORM 

TO EOUICONTROL SPACE ANO COMPUTE EIGENVALUES OF OVD 

0E60271O 

C 



0EG0272O 



CALL VMULFF(0,V.N,N.N. 10. 10.0V. 10.IER) 

0EQ0273O 



CALL VMULFF(DV.O.N.N.N. 10. 10.0VO. 10.IER) 

0EG0274O 



CALL VCVTFS(0VO.N, 10.0V0SYM) 

0EQ0375O 


ooooo non 
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FILE: OEQCON FORTRAN A 


VM/SF CONVERSATIONAL MONITOR SYSTEM 


CALL ElCRS(0VOSVM.N,I0O8.EV.Z. 10.WK.IER) 
IF((ABS(EV(1)).LT.T0L).0R.(EV(1).LT,0.)) 00 TO T6 

COMPUTE DEGREE OF CONTROLLABILITY 

VS-S0RT(EV(1)**N) 

PROOEV-t.O 
00 706 1-1. N 

PROOEV-PROOEV* EV< I } 

706 CONTINUE 

VE-SORT(PROOEV) 

POWER- 1.0/FL0AT(N) 

OEGCON- ( VS«( VS/VE ) • ( VE - VS ) ) —POWER 
GO TO 80 

76 OEGCON-0. 

• STORE OC IN APPROPRIATE ARRAY; SEARCH FOR MAXIMUM DC 

» AND RECORD ITS LOCATION. MAGNITUDE AND MAXIMUM ANO 

* MINIMUM CONTROLLABLE DIRECTIONS 

80 IF(NA.NE.4) GO TO 83 

0C4(IACT4. IACT3. IACT2. lACT 1 )-OEGCON 
IF (OEGCON. GT.OCMAX) GO TO 80S 
GO TO ISt 

80S lF((IACT1.E0.1ACr2).0R.(IACT1.E0.IACT3).0R. 

& ( I ACT 1 . EO . I AC T4 ) . OR . ( I ACTS . EO . I ACT3 ) . OR . 

& (IACT2.EO.:^CT4).OR.(IACT3.EO.IACT4)) GO TO 1S1 

OCMAX -OEGCON 
LMAX-LOC 
00 807 I-1.N 

UMIN( I )-Z( I . 1 )/0( I . I ) -SORT(EV( 1 ) ) 

UMAX( I )-Z( I .N)/0( I . I )*SORT( EV(N) ) 

807 CONTINUE 

GO TO 1S1 

83 IF(NA.NE.3) GO TO 87 

0C3( 1ACT3, 1ACT2. 1 ACT 1 ) -OEGCON 
IF (OEGCON. GT. OCMAX) GO TO 83S 
GO TO 1S1 

83S IF((IACT1.EO.IACT2).OR.(IACT1.E0.1ACT3).OR. 

8 (1ACT2.E0.IACT3)) GO TO 1S1 

OCMAX -OEGCON 
LMAX-LOC 
00 847 1-1. N 

UMIN( I )-Z( I . 1 )/0( I . I )-SORT(EV( 1 ) ) 

UMAX( I )-2( I .N)/0( I . 1 )*SORT( EV(N) ) 

847 CONTINUE 

GO TO 1S1 

87 IF(NA.NE.2) GO TO 89 

0C2( 1ACT2. IACT1 }-0EGC0N 
I F( OEGCON. QT. OCMAX) GO TO 87S 
GO TO tSI 

87S IF( IACT1.E0.IACT2) GO TO 1S1 

OCMAX -OEGCON 
LMAX-LOC 
00 877 I-1.N 


DEQ02760 

0EGO277O 

DEG02780 

0EOO279O 

0EQO28OO 

0EGO2B1O 

0EGO282O 

0EQO283O 

0EG0284O 

0EGO28SO 

0EGO286O 

0EG0287O 

0EG02880 

0EG02890 

DEGU2900 

0EGO291O 

0EGO292O 

0EG02930 

0EGO294O 

0EGO29SO 

0EGO296O 

0EGO297O 

0EGO2980 

0EQO299O 

0EG03000 

0EGO3O1O 

0EG03020 

0EGO3O3O 

0EGO3O4O 

0EGO3OSO 

OEG03060 

0EGO3O70 

0EGO3O8O 

0EQO3O9O 

0E6O31(X> 

0EGO31 10 

0EGO312O 

0EGO313O 

0EGO3140 

0EGO315O 

DEG03160 

0EGO317O 

0EGO318O 

0EGO319O 

0EGO32(X) 

0EG03210 

0E0O3220 

0EQO323O 

0EQ0324O 

0EG032SO 

DE003260 

0EG03270 

0EGO3280 

0EQO329O 

0EOO33OO 


\ 



o o o 


FILf: OEGCON FORTRAN A 


VM/SP CONVCRSATIONAL MONITOR SVSTEM 


UMIN( I )-Z( I . 1 )/D( 1 . I )*SORT(EV( 1 ) ) 

UMAX( I )-2( 1 .N)/0( I . I >*SORT(EV(N) ) 

877 CONTINUE 

00 TO 181 

89 OCKUCTD-OEOCON 

IP ( OEGCON. GT.DCMAX) GO TO 898 
GO TO 151 

895 OCMAX-OEGCON 

* LMAX-COC 

00 897 1*1, N 

UMIN( I )*2( I . 1 )/0( I . I )*SORT( EV( 1 ) ) 
UMAX(I)-Z(I.N)/0(I.n*SORT(EV(N)) 

897 CONTINUE 

181 CONTINUE 

IP(NA.LT.2) GO TO 203 
161 CONTINUE 

IF(NA.LT.a) GO TO 203 
171 CONTINUE 

IF(NA.LT.4) GO TO 203 
181 CONTINUE 

• COMPUTE average DC AND SEARCH FOR MAXIMUM 

203 OCMAXF-0. 

IF(NA.NE.4) GO TO 300 
00 250 I*ISTART.1FC00E 
00 250 O-ISTART.IFCOOE 
00 250 K-ISTART.IFCOOE 
00 280 L-ISTART.IFCOOE 

CALL PM4EXP (OCA , I FCOOE . TAU . TOP . I . J . K . L . OCA VE ) 

LOC-I • 10**6*J* 10**4aK* 100+L 
WRITE! 10.228) L0C.0C4(1, J.K.D.OCAVE 

FORMAT! ‘LOCATION- 18. SX. ‘OC- ‘ , E 1 1 . 4 . 5X . 'OCAVE- ' . E 1 1 . 4 ) 
IF!OCAVE.GT.OCMAXF) 60 TO 230 
GO TO 250 

IF!!IEO.J).OR.(I.EO.K).OR.!l.EO.L).OR. 

!O.EO.K).OR.!J.EO.L).OR.(K.EO.L)) 60 TO 250 
LMAXF ■ I • 10- -e^d- 10- -A+K- 10--2*L 
OCMAXF-OCAVE 

250 CONTINUE 
GO TO 700 

300 IF!NA.NE.3) GO TO 400 
00 350 I-ISTART.IFCOOE 
00 350 d- 1 START. I FCOOE 
00 350 K-ISTART.IFCOOE 

CALL PM3EXP ! 0C3 . I FCOOE . TAU . TOP . I . d . K . OCAVE ) 

IFIOCAVE dT.OCMAXF) 60 TO 330 
GO TO 380 

330 IF(!l.EO.d).OR.(I.EO.K).OR.!d.EO.K)) GO TO 380 

LMAXF -I FCOOE -10- -6- 1* 10* 'A^d* 100'^K 
OCMAXF-OCAVE 
380 CONTINUE 
GO TO 700 

400 IF!NA.NE.2) GO TO 800 
00 450 1 -I ST ART. I FCOOE 


220 

225 

230 

& 


0EGO331O 

0EG03320 

OEG03330 

OEQ03340 

0E003380 

OEQ03360 

OEQ03370 

0EQO338O 

DEQ03390 

0EQ03400 

0E0034 10 

0EGO3420 

0EGO343O 

0EGO344O 

0EGO348O 

0EG03460 

0EGO347O 

0EQO348O 

0EGO349O 

0EGO3SOO 

0EGO38 10 

0EGO382O 

0EGO3B3O 

0EGO354O 

0EGO3S8O 

0EGO356O 

0EGO357O 

0EG03580 

0EQO3S9O 

0EGO36OO 

DE0O361O 

0EGO362O 

0EQO363O 

0EGO364O 

0EGO368O 

0EGO366O 

0EGO367O 

0EQ0368O 

0EGO369O 

OEG03700 

0EQO371O 

0EGO372O 

0E603730 

OEQ03740 

0EGO375O 

0EG03760 

0EQO3770 

0EGO378O 

0EGO379O 

0EQO3800 

OEG03810 

OE 0038 20 

0EGO383O 

0EG03840 

0EGO385O 


O O O U U U V 
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FILE; OEGCON FORTRAN A 


VM/SR CCmVERSATIONAU MONITOR SYSTEM 


WRITE! 10.433) L0C,0C3(I.J),0CAVE 

433 FORMAT! 'LOCATION-' . It.SX. '0C> ' ,E t1 .4.8X. 'OCAVE-' , E 11 . 4) 

IF!OCAVE.GT.OCMAXF) GO TO 430 
GO TO 4S0 

430 IFU.EO.J) GO TO 4S0 

LMAXF - 1 FCOOE • 10- *6-t> I FCOOE • 10- -44 1 • 100*0 
OCMAXF-OCAVE 

ABO CONTINUE 
GO TO 700 

500 00 550 1- 1ST ART. I FCOOE 

LaC-IFCOOE-lO-«6*tFCOOE»10--4*IFCOOE- 100*1 
WHITE! 10.535) LOC.OCi!!) 

525 FORMAT! 'LOCATION-' .18. 5X. '0C-' .E11.4) 

550 CONTINUE 

* OUTPUT DC'S. LOCATIONS. AND PRINCIPAL DIRECTIONS 

700 IFlNA.NE.4) GO TO 711 

WRITE ! 10 . 70S ) DCMAX . LMAX .DCMAXF . LMAXF . 
a !UMIN!1).1-1. 10).(UMAX!I).1-1. 10) 

705 FORMAT! IX. 'MAX OC FOR 4 OPERATIONAL ACTUATORS IS'. El 1.4/ 
a 'ANO THE LOCATION IS '.18// 

a 'MAX OC FOR 4 FAILING ACTUATORS IS'. El 1.4/ 

a 'ANO THE LOCATION IS '.18// 

a 'UMIN-'/5!E11.4.5X)/5!E11.4.5X)// 

a 'UMAX-'/5!E11.4.SX)/S(E11.4.SX)//) 

GO TO 1000 

711 IF!NA.NE.3) GO TO 714 

WRITE!IQ.71S) OCMAX , LMAX . DCMAXF . LMAXF . 

a !umin!i).i-i. io). !umax(i).i-i. io) 

715 FORMAT! IX. 'MAX OC FOR 3 OPERATIONAL ACTUATORS IS'. El 1.4/ 
a ‘ANO THE LOCATION IS '.18/ 

a 'MAX OC FOR 3 FAILING ACTUATORS IS'. El 1.4/ 

a 'ANO THE LOCATION IS '.18// 

a 'UMIN-'/5!E11.4.5X)/5(E11.4.5X)// 

a 'UMAX-'/5!E11.4.5X}/S!E11.4.SX)//) 

?:C TO 10O0 

714 IF(NA.NE.2) GO TO 731 

WRITE! 10.720) LMAX. DCMAX. LMAXF. OCMAXF. 

4 !UMIN! 1 ) . I- 1 . 10) . !UMAX! I) . 1- 1 . 10) 

720 FORMAT! //' LMAX- '. 18. 10X. 'OCMAX- ' . E 1 1 . 4/ 

a 'LMAXF-' .18. 10X. 'OCMAXF-' .Eli. 4// 

a 'UMIN-'/5!E11.4.5X)/S!E11.4.5X)// 

a 'UMAX-'/5!E11.4,5X)/5!E11.4.5X)//) 

GO TO 1000 

721 WRITE!I0.73O) LMAX. OCMAX. 

a (umin( 1 ).i-i.io).!umax(i).i-i.io) 

730 FORMA' .MAX-' . 18. 10X. 'DCMAX-' .El 1.4// 
a 'UMIN-'/5!E11.4.5X)/5!E11.4.SX)// 

a 'UMAX-'/S!E11.4,SX)/5!E11.4.SX)//) 

> PLOT OF OC VS. ACTUATOR POSITION FOR 

> 1 FIXED ANO 1 VARIABLE ACTUATOR 

1000 00 1002 1-1.21 


0EQ0386O 
OEQ03870 
OEQ03880 
0EG03890 
0EG0390O 
0S0O391O 
OE003820 
DEQ0383O 
0EQ03840 
0EQ0385O 
0EQ0398O 
0EG0397O 
0EGO398O 
0EG0399O 
OEQ04000 
0EGO401O 
0EQO4O2O 
0EGO403O 
DE00404O 
0EQO409O 
0EQ040SO 
0EQO407O 
OE004080 
0EQ0409O 
0EGO4 100 
0EQO411O 
0EQ0412O 
0EQO413O 
0EQO414O 
0E00419O 
0EQ0416O 
DE00417O 
OEQ04180 
0EQ0419O 
0EQO42OO 
0EQO421O 
0EQO433O 
0EQO423O 
0EQO434O 
0EQO42SO 
DEQ04260 
0EQO427O 
0EQO428O 
0EQO429O 
0EQ04300 
0EQO431O 
0EG0432O 
0E0O433O 
0EQ04340 
0EGO439O 
0EQC436O 
0EQO437O 
0EQ0438O 
0EQO43BO 
0EGO44OO 



PILE: OEOCON EORTRIN A 


VM/SR CONVERSATIONAL MONITOR SYSTEM 


XARRAV\ I )»RL«EL0AT( I • 1 )/20. O 
1002 CONTINUE 

00 1004 I ■ 1 . 1 1 

YARRAY(I)-0C2(IPIX.I) 

1004 CONTINUE 

00 1009 1-12.21 

VARRAY(I)-0C2(IEIX.22-I) 

1009 continue 

CALL SCALC(XARRAY.«.0,2M) 

CALL SCALE(YARRAY.4.0.21. 1 ) 

CALL AXI$(0. .0. .‘SEAM R0SIT10N( FT ) ' . - 17. 6 .0.0.0. 

S XARRAV(22),XARRAY(23)) 

CALL AXIS(0. .0. .‘OEOREE OF CONTROLLABILITY' .♦29.4,0.90.0. 

S YARRAY(22).Y4RRAY(23)) 

CALL LINE(XARRAY. y ARRAY. 21 . l.ol.S) 

CALL SYMB0L(0. 9. 9. 0.0.21. 'DEGREE OF CONTROLLABILITY' .0.0. 29) 
CALL SYMBOL! 1.0.4,9.0.21. 'FOR A FREE*FREE BEAM' .0.0.20) 

CALL ENOPLT( 12.0.0.0.999) 

STOP 

2000 WRITE! 10. 100) 

100 FORMAT! 'THE MINIMUM E-VALUE IS ZERO') 

00 2001 1 - 1.10 

WRITE! 10. 2002) EV!l) 

2002 FORMAT! 'EV-'. Eli. 4) 

2001 CONTINUE 
r.TOP 
END 

C 

C- MODAL AMPLITUDE AT X FOR SIW>LY-SUPPORTEO DEAM 
C 

REAL FUNCTION PH12!X.M00E.BM.BL) 

DATA PI/3. 141992994/ 

PHI2-S0RT! 2 .0/BM)-SIN! FLOAT! MODE ) -PI -X/BL ) 

RETURN 

END 

C 

C« MATRIX ! ARRAY) TjMES VECTOR (V) 

C 

SUBROUTINE MATVEC! M.N. ARRAY . V . RET ) 

DIMENSION ARRAY!M.N).V!N) RET!M) 

DO 10 1-1. M 
RET!l)*0. 

00 10 J-1.N 

10 ret!i)-pct':>*array!i. j)-v!j) 

RETURN 

END 

C 

C* ADDS MATRIX 8 TO A 
C 

SUBROUTINE HATA00!N. A.B.RET) 

DIMENSION A!N.N).B!N.N).RET!N.N) 

00 10 I-1.N 
00 10 J-1.N 

10 RET!I.U)-A(I.J)yB!l. J) 

RETURN 


0CaO441O 

0EQO4420 

OE004430 

OE004440 

0EQO44BO 

OE004490 

0C004470 

0E0O44BO 

0EQO449O 

0E'JO490O 

0EQO491O 

0EQO492O 

0EO0493O 

0EB04940 

DEQ04990 

DE004S9^> 

0EQO4970 

OE004980 

0EQO499O 

OEQ04900 

0EQ04910 

OE004920 

DE004630 

DEQ04B40 

0EQ0495O 

0EQ0496O 

0EQ04970 

0E0O46S0 

0EQ0469O 

0EGO4VOO 

0EOO471O 

0EQO4720 

0EQO473O 

0EGO4740 

0E004790 

0EQO476O 

0£004V70 

0EOO478O 

DEQ047B0 

0EOO4BO0 

OE004910 

DE004B20 

DEa04830 

01.004840 

0EOO4B9O 

OEG04860 

0E004S70 

OEQ04880 

0EQO489O 

OE004900 

0EGO491O 

OE004920 

0EQO493O 

0EQO494O 

DEG04990 


UUU uuu uuuuu uuu 


DEOCON POKTKAM A 


VM/S^ COfViftSAriONAL MONITOR SYSTEM 


ENO 

• SUBTRACTS MATRIX B FROM A 

SUBROUTINE MATSUB(N. A.B . RET ) 

DIMENSION A(N.N). BIN. N). RETIN, N) 

00 10 1-t.N 
DC 10 0-1.N 

10 RET(l,J)-A(I.U)-B(I.d) 

RETURN 

END 

• MOOAL AMPUTUOE AT X FOR FREE-FREE BEAM 

REAL FUNCTION PHI ( X .BETA ,BL ) 

ALP-BETA<BL 

SH-0.3*IEXP(ALP)-EXP(-ALP)) 

CH-0. 5«( EXP( ALP )^EXP( -ALP ) ) 

A-(SH«SXN(ALP))/(CH-COS(ALP>) 

PHI>0.S*(EXP(BETA«X)«EXP(-BETA*X))^C0S(BETA*X)- 
A A*(0.B*(EXP(BETA*X)-EXP(-BETA*X))*S1N(8ETA.X)) 

RETURN 

ENO 

• COMPUTES AVERAGE EXPECTED PERFORMANCE MEASURE FOR 

• 4 COMPONENTS ASSUMING EACH HAS SAME EXPONENTIAL 

• DISTRIBUTION OF TIME TO FAILURE 

SUBROUT I NE PM4 EXP( PW. I FCOOE . TAU . TOP . L 1 . L3 . L3 . L4 . PMAVE ) 
DIMENSION PM( 13. ia. 13. 13) 

tt-top/tau 

PT1-< 1.0/(4.0«TT))*( 1.0-BXP(-4.0*TT)) 

PT33-1,0/( 13.0*TT)-( 1.0/( 13 .0*TT ) )*(4.0-3.0*EXP( -TT ) ) • 

A EXP(-3.0«TT) 

PT611-1.0/( ia.O»TT)-( 1,0/( 13.0*TT))*(6.0-B.0*EXP( -TT)^ 

A 3.0»£XP(-2.0*TT))»EXP(-3.0»TT) 

PT 13 18-1 .0/(4.0*TT)-( 1 .0/(4 ,0*TT) )*(4 .0-6 ,0-EXPI -TT)^ 

A 4.0*EXP(-a.0-TT)-EXP(-3.0-TT))*EXP(-TT) 

PT16-1.0-( 1,0/( ia.0*TT))-(-38.O'4B.O-EXP(-rT)+ 

A 36 . 0*EXP( -3 . 0-TT ) - 16 . 0*EXP( -3 . 0*TT )* 3 . 0*EXP( -4 . 0*TT ) ) 

PMAVE -PT 1 -PM( L 1 . L3 . L3 . L4 )>PT3S-( PM( IFCOOE . L3 . L3 , L4 )♦ 

A PM( L 1 . IFCOOE . L3 . L4 ) ♦PMl L 1 . L3 . I FCOOE . L4 )♦ 

A PM( L1,L3.L3. I FCOOE ))^PT611*(PM( IFCOOE. IFCOOE. L3,L4)4> 

A PM( I FCOOE . L3 . I FCOOE . L4 ) ♦PM( I TCOOE . t3 . L3 . I FCOOE ) ♦ 

A PMI L t . I FCOOE . I FCOOE , L4 ) ♦PM( L 1 . 1 FCOOE . L3 . 1 FCOOE ) ♦ 

A PMI L t. L3. IFCOOE. I FCOOE ))«PT 13 IS- (PM( I FCOOE. IFCOOE. 

A I FCOOE. L4) ♦PMI t .'CODE . IFCOOE . L3 . IFCOOE )♦ 

A PMI I FCOOE . L3 . 1 FCOOE . 1 FCOOE PMI L 1 . 1 FCOOE . I FCOOE . IFCOOE ) ) 

RETURN 
ENO 

SAME AS PM4EXP EXCEPT FOR 3 COMPONENTS 

SUBROUTINE PM3EXP( PM . I FCOOE . TAU . TOP . L 1 . L3 . PMAVE ) 

DIMENSION PMI 13. 13) 


0EQO496O 
0E004B70 
0EQO4fBO 
0EQO46BO 
OEQ09000 
0EQOB01O 
OE008030 
0EGO903O 
OEGOSORO 
OEOOSOSO 
0EG0S06O 
0EOOSO7O 
OEGOSOBO 
OEOOSOBO 
OEOOBtOO 
0E00911O 
0EQO813O 
DE00S130 
DEQ0S14O 
OEOOB1SO 
0EQ0S16O 
DE00S17O 
OEQOStBO 
0EQ09190 
0EQ09300 
0E00831O 
DEQ08330 
0E00B33O 
0E00S34O 
OEOOS3SO 
0CQ0S36O 
0E00937O 
0EG0836O 
0E0OS39O 
OE 008300 
0EOOS31O 
0EGO833O 
OEOOS330 
DEQ0934O 
DEQOS3SO 
0EQO936O 
0EQOS37O 
0EGOS38O 
0EG0839O 
0EGOS40O 
OEOOS4 10 
0EGOS43O 
0EGOS43O 
0EGOS44O 
DEGOSASO 
0EGO546O 
0EGOS47O 
DE 0084 BO 
OEQOS490 
OEGOSSOO 
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I* I LI: OCOCON POKTRAN A 


VM/SR convcrsational monitor SVSTCM 


TT-TOR/TAU 

RTt«( 1.0/(a.0*TTn«(1.0-fXR<-a.0*TT)) 

RTaa«n.O/TT)«(( 1.0-IXR(-TT))-0.8*( I.O-IXR(-a.O»TT))) 

RMAVC -RT 1 »RM( L 1 . L3 )*RTaa* ( RM( 1 PCOOl . Li ) *RM< L 1 . 1 RCOOf ) ) 

RETURN 

END 

C 

C* SAME AS RM4EXR EXCBRT FOR 3 COMRONENTS 

C 

SU8RCUT I NE RM3EXR( PM . I FC(X>E , TAU . TOR . L 1 . L3 . L3 . PMA VE ) 

DIMENSION RM(ia.ia.ia) 

TT-TOR/TAU 

RTl-( 1.0/(3.0*TT))«( 1.0-EXR(-3.0*TT)) 
RTa4*(1.0/(a.0*TT))»(1.0-EXR(-a.O*TT))-( 1.0/0. 0-TT))* 

A (1.0-EXP(-3.0*TT)> 

RT57*( 1.0/TT)»( t.O-EXR(-TT))>( 1.0/TT)*{ 1.0-EXR(-a.0*TT))^ 

A (1.0'O.0*TT)i*( 1.0-EXP(-a.0«TT)) 

RMAVE -PT 1 •RM( L 1 . La . L3 ) 

A ♦PTaA* (PM( I FCOOE . L2 , 1 3 )-«>PM( L 1 . 1 FCOOE . L3 )^RM( L 1 . L2 . IFCOOE ) ) 

A ♦PT87«>(PM( IFCOOE . IFCOOE . L3 )>PM( IFCOOE . L2 . IFCOOE )♦ 

A PM(L 1.1 FCOOE. I FCOOE)) 

RETURN 

END 

C 

C* COMPUTES DIAGONAL SOLUTION BLOCKS 0» M (TYPE III) 

C 

SUBROUTINE 0IAO(0T.SIO1 .0M1 .01 1 .012.022 .AM1 1 , AM12. AM22 ) 

DATA EPS/0.000001/ 

IF{AeS(SIG1).LT.ERS) GO TO S 
A-(01 1^J22)/(4.0«SIQ1) 

B-(O.S*Siai*D12-O.2B*OM1*(n22-011))/(0M1*0M14SIQ1*SI(;i) 
C-(O.S*0Ml*012^O.2S*SIG1*(D22-01 1))/(0M1*0M14SIQ1«SIQ1 ) 
S--2.O»S101»0T 
AR0--2.0*0M1*0T 

AM11-A*( 1.0*EXR(S))*B«EXP(S)»SIN(ARG)-C«( 1.0-EXR(S)*C0S(AR0)) 
AM12->C*EXP(S)*SIN(AR0)4-8*( 1.0-EXR(S)*C0S(AR0)) 

AM22-A*( 1.0-EXR<S))>B*EXP(S)*SIN(ARG)^C*( 1 -0-EXP(S)*COS(ARO) ) 

GO TO 10 

S A>0.5*(022«01 1 ) 

B-(022-011)/(4.0*0M1) 

C-012/(2.0*OM1 ) 

ARa-'2.O*0M1*0T 

AMI 1--A*(-0T)^B*SIN(ARG)-C*( I.O-COS(ARG)) 
AM12--C*SIN(ARG)-8*(I.O-COS(ARG)) 

AM22-'A*(-0T)-B*StN(ARG)^C*( 1 .0-C0S( ARG) ) 

10 RETURN 
END 
C 

C« COMPUTES OFF-OIAOONAL SOLUl I(M BLOCKS OF ^ (TYPE IV) 

C 

SUBROUTINE OFOIAQ(DT.SIG1.SIG2.0M1,OM2.011.D12. 

A 021.022.AM11.AM12.AM21.AM22) 

SIGT-SIGUSIOa 
S-SIGT*( -OT) 

ARaPx(0M2«0M1)«(-0T) 


OBQ08B10 
OE 008820 
OE008530 
0EQ0884O 
0.308880 
0EQ088A0 
0EGO8870 
0E0O88AO 
0EQOSS8O 
OEOO^AOO 
OEOOSA 10 
0E0OS820 
0EGO8C3O 
0E0O8A4O 
OEOOSASO 
0EQOSB80 
OE008A70 
OEGOSAAO 
OEQOS8BO 
OEG08700 
0EGO871O 
0EQO872O 
0EOO873O 
OEGOS74C 
DE0O87S0 
0EOO87BO 
0CQOS77O 
0E0OS780 
0EQOS79O 
0EQOS8O0 
OEOOSA 10 
01008820 
0E00S830 
OE008840 
0EQ08880 
OEQOSeSO 
0E0O987O 
0E0088AO 
0E00S890 
0EGCS900 
OEGOS810 
OEOOS920 
0EG0893O 
0EQ09940 
0E0059SO 
0E0O99AO 
0EGO997O 
0EGO99AO 
OEQ08990 
OEOOeOOO 
OEG06C10 
OE 006020 
0S006030 
OEG06040 
0EQO6O9O 


3-12 


I 

FILt: OEQCON EOKTRAN A 


VM/SP CONVERSAriON/L MONITOR SYSTEM 


ARQM- ( OM3 -OR I ) • ( -OT ) 

A--((OMa'OM < )*(01t^022)^SIQT*(Ol2‘Oai))/ 

« (3.0'>((CM2-OMt)**2^SiaT**2)) 

■■((0M2*QM«)«(0at*012)«SiaT«(0n4O22))/ 

» (2.0*((OM2'OM1)*«24^SIQT**2)) 

C-*((OM2«OMI ) *(011-022 >^SiaT«(D2UO 12))/ 

» (2.0*((0M2«0M1 )**2>SIOT**2)) 

0--((OM2«OM1 )• (021^0 12 )^SIQT*(022-011 })/ 

S (a.O*^(OM2tOM1)**2^SiaT*-2)) 
AM11>l«0«EXP(S)*(A*SIN(AR<2M)-a*COS(ARQK)^ 

• C*SIN(ARQP)-0*COS(Aroi»)) 
AM12--A-C*EXP(S)*(R*SIN(AR0N)*A«C0S(ARGM)4- 

• 0*S1NURQP)«C*C0S(XRQP)) 
AMai-A-C>EXP(S)«(-B*SIN(AINM)-A*COS(ARaM)^ 

« 0*SIN(ARQP)4^C*C0S(ARQP)) 

AM22-R-0«EXP(S)*(A*SIN(ARaM)-B*C0S(ARQM)- 

• C*SIN(AROP)«0*COS(AROP); 
RETURN 

END 


OEOOB060 

0E00B070 

OEOOBOBO 

OEOOBOBO 

0E0OB1OO 

0EQ0B110 

0E00B120 

0E00B130 

0EO0B14O 

OEOOB150 

OEGOBieO 

0EGOBI70 

OEQOBtiO 

0E00B190 

0E0OB20O 

OEOOBaiO 

OEOOB220 

0EG0B230 

0EOO«2«O 


V 


